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Abstract — The rapid proliferation of wireless systems malies 
interference management more and more important. This paper 
presents a novel cognitive coexistence framework, which enables 
an infrastructure system to reduce interference to ad-hoc or peer- 
to-peer communication links in close proximity. Motivated by the 
superior resources of the infrastructure system, we study how its 
centralized resource allocation can accommodate the ad-hoc links 
based on sensing and predicting their interference patterns. 

Based on an ON/OFF continuous-time Markov chain model, 
the optimal allocation of power and transmission time is formu- 
lated as a convex optimization problem and closed-form solutions 
are derived. The optimal scheduling is extended to the case where 
the infrastructure channel is random and rate constraints need 
only be met in the long-term average. Finally, the multi-terminal 
case is addressed and the problem of optimal sub-channel 
allocation discussed. Numerical performance analysis illustrates 
that utilizing the superior flexibility of the infrastructure links 
can effectively mitigate interference. 

Index Terms — Cognitive Radio; Resource Allocation and Inter- 
ference Management; Standards Coexistence; Dynamic Spectrum 
Access; 

I. INTRODUCTION 

The rapid growth of wireless networks makes interference 
an important performance impediment and motivates a careful 
study of coexistence. In unlicensed bands, where there is a lack 
of coordination among heterogeneous technologies, traditional 
forms of coexistence are based on statically separating sys- 
tems, for example by having them operate in disjoint frequency 
bands. However, such static approaches are inefficient and not 
able to accommodate the projected growth in deployments. 

The high interference conditions which are prevalent in 
unlicensed bands are in contrast to a fairly low average 
utilization, which results from the burstiness of traffic that is 
supported by these systems and the random medium access 
behavior typically employed. This contrast motivates dynamic 
approaches to interference management which can leverage 
unused spectrum opportunities. 

Cognitive radio presents a new framework for analyzing 
this problem. Reconfigurable radio platforms enable a dynamic 
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adaptation of transmission parameters based on monitoring 
the radio environment. In this way it is possible to reduce 
interference by detecting and predicting temporal activity 
patterns and avoiding transmission overlaps. Ultimately this 
may lead to a new cognitive coexistence paradigm. 

This paper addresses a special case of cognitive coexistence 
involving two different types of networks: an infrastructure 
(IS) wide area network that shares spectrum with local, ad- 
hoc (AH) or peer-to-peer systems. Motivated by the superior 
communication resources of the IS system, we analyze how its 
flexible, centralized resource allocation can accommodate the 
AH links based on sensing and predicting their interference 
patterns. Despite adapting its resource allocation based on 
sensing results, the IS system minimizes interference to the 
AH network subject to maintaining a specified quality-of- 
service (QoS) level for its users. 

This approach is different from typical dynamic spectrum 
access (DSA) problem formulations in which a secondary 
system exploits spectrum opportunities left over by a primary 
system, subject to the constraint that no significant interference 
is created (see HI for a review). While both approaches are 
hierarchical, spectrum property rights in DSA mandate that 
secondary users adapt to spectrum licensees. In contrast, in 
cognitive coexistence it is possible to use the flexibility of the 
primary system to accommodate a secondary (perhaps lower 
priority) system while maintaining a desired primary system 
performance level. This could be viewed as a "best-effort" 
approach toward interference management and coexistence. 

This framework is relevant to a number of practical se- 
tups. For example, the convergence of wide and local area 
networks has received increasing interest. The coexistence 
of IEEE 802.16 and IEEE 802.11 systems in unlicensed 
bands is such an example of practical importance 121, Q. 
Similarly, incorporating peer-to-peer connectivity into cellular 
networks [4| or accommodating femto-cell base stations Q 
are potential scenarios where this work could be applied. Yet 
other applications arise in the military domain, where the 
coexistence of high and low priority finks is a fundamental 
concern JS). 

A. Main contribution 

This paper addresses the coexistence of two wireless net- 
works that operate on different spatial scales: a longer range 
IS network that interferes with local AH links. We address 
the question of whether the IS network can accommodate sur- 
rounding AH links by allocating power and transmission time 
judiciously. Specifically, we make the following contributions: 
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(b) Time/frequency setup 

Fig. 1: System setup. An infrastructure link allocates power and transmission time such as to minimize interference to close-by 
ad-hoc networks. The interference- aware resource allocation is based on detecting and predicting the ad-hoc system's temporal 
activity. 



• The problem of optimal power and transmission time 
allocation is formulated as a convex program and the op- 
timal frame-level solution is derived. Based on optimality 
conditions, a solution algorithm with guaranteed conver- 
gence and low complexity is introduced. The structure of 
the optimal solution is studied to provide further intuition. 

• The frame-level problem is relaxed to the case of average 
rate constraints, in which statistical knowledge of the 
activity patterns of the AH users and the IS channel co- 
efficients are used to allocate resources in both frequency 
and time. 

• Finally, a scenario with multiple IS users is considered 
in which the IS base station allocates sub-channels based 
on average interference metrics and users perform op- 
timal power and transmission time allocation based on 
sensing results. A comparison with conventional sub- 
channel allocation methods shows that heuristics may 
yield performance close to optimal. 

These contributions are corroborated by numerical results 
which demonstrate that judicious transmission time allocation 
can mitigate interference effectively. 

B. Related work 

Optimal resource management in multiuser multicarrier 
wireless systems has been well-studied for both downlink and 
uplink cases; see Q, lEI, lH for an overview of the topic. 
In cognitive radio networks, optimal resource allocation is 
more challenging because it needs to incorporate interference 
constraints, which protect the primary system from harmful 
interference. Typical formulations aim at finding a power and 
sub-channel allocation which maximizes the throughput of the 
cognitive radio system while meeting interference and power 
constraints. Recent work in this area includes ifTOl . ifTTI . In 
addition to meeting interference constraints, spectrum sharing 
and self-coexistence within the cognitive radio network also 
need to be addressed. Contributions in this area include lfT2l . 
ifTsl . lfT4l . ifTsl . Within the framework of cognitive coexis- 
tence, the optimal power allocation based on knowledge of the 



interference channel has been addressed in lfT6l by the authors 
of this paper. Previous work on improving the coexistence 
among local and personal area networks includes ifTTl . ifTSll 
in which a cognitive frequency hopping protocol is derived 
based on temporal activity models. To the best of the authors' 
knowledge, interference-aware power and transmission time 
allocation based on predicting temporal activity patterns has 
not been addressed before. 

C. Organization and notation 

The rest of this paper is organized as follows. After intro- 
ducing the problem setup in Sec. |ll] the optimal frame-level 
allocation is derived in Sec. |lll] The results are extended to 
the average rate case in Sec.|lV]and sub-channel assignment is 
considered in Sec. |Vl Throughout this paper notation is fairly 
standard. Vectors are typeset in boldface. For an event X, the 
indicator function l^x] is equal to one if X occurs and zero 
otherwise. The notation (•)+ is used to abbreviate max{0, ■}. 

II. Problem Formulation 
A. System setup 

The system setup is shown in Fig. [T| We consider an IS 
system, which consists of a base station and a single client (the 
multi-terminal case will be addressed in Sec. |V]i. The uplink 
transmissions of this client may strongly interfere with local 
transmissions of one or multiple AH networks surrounding the 
client. For this setup, the problem of optimally assigning power 
and transmission time at the IS client, such as to minimize 
interference to the AH links is analyzed. The time/frequency 
behavior of both systems is shown in Fig. [lb] 

Ad-hoc network: The AH network consists of a set of AH 
nodes which operate in a frequency band that overlaps with the 
IS system. As depicted in Fig. [lb] there can be multiple AH 
networks which operate in non-overlapping bands that each 
overlap with a certain set of IS sub-channels. It is assumed that 
the partitioning of the AH bands is fixed and that the temporal 
activity of different bands is statistically independent. 
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We model the time behavior of each AH band by a two-state 
ON/OFF continuous time Markov chain (CTMC). The holding 
times in both ON and OFF state are exponentially distributed 
with parameters ^ for the ON state and A for the OFF state. 
Therefore, if an AH link is detected to be in a certain state at 
time to, then its transition matrix for time to + t is given by 



P(7 



1 



(1) 



which follows directly from the definition of a CTMC |fT9l 
p. 391]. Therefore, the probability of an AH link being ON 
at time to + r, conditioned on having it observed in the ON 
(OFF) state at time to is given by the lower right (upper right) 
entry in the matrix above. 

Modeling AH links based on a two-state CTMC approx- 
imates the carrier sense random medium access typically 
employed in such systems. This modeling approach has been 
used in related publications ifTSl . Il20ll and solidified by 
a measurement-based analysis of WLAN traffic ifTTl . The 
CTMC assumption strikes a good tradeoff between model 
accuracy and the analytical tractability that is needed in the 
subsequent sections. 

Infrastructure system: The IS system operates in the same 
frequency band as the AH network and evolves in frames of 
fixed duration T. At the beginning of each frame, spectrum 
sensing is used to detect the ON/OFF activity of the AH bands, 
and based on the sensing outcome, power and transmission 
time are assigned; see Fig. [Tbl 

Throughout the paper we assume perfect sensing, that 
is, sensing outcomes are always accurate and the overhead 
associated with sensing is neglig ibleQ Due to the proximity 
of IS and AH systems, the detection task is conceptually 
similar to the carrier sensing employed in systems such as 
IEEE 802.11. 

Based on the sensing result at the beginning of each frame, 
the IS system allocates power and transmission time on a sub- 
channel basis. This is conceptually similar to the allocation of 
time/frequency resource blocks in broadband cellular systems 
based on OFDMA. The case where a subset or even all sub- 
channels need to share the same timing allocation (for example 
when transmissions in the entire band can only be turned on or 
off) has worse performance in general. Nevertheless, a similar 
solution approach remains applicable. 

The IS system minimizes interference subject to maintaining 
rate requirements for its client. The rate that is supported by a 
specific sub-channel is modeled based on a channel capacity 
formulation, 

Pn log 1 + n — = > p„ log 1 + 

V PnNo ) ^ \ Pn 

(2) 

where p = [pi, . . . ,pjv]^ denotes the power allocation, 
p = [pi, . . . , pn]'^ represents the transmission time allocation, 

'in contrast to many DSA setups, where very weak signals need to be 
detected, the proximity of AH terminals to the IS client leads to moderate 
to high average signal-to-noise-ratios. This facilitates the sensing task and 
enables us to employ simple methods such as energy detection. 



A^o is the noise power, k a normalization factor, and /3„ is 
introduced for notational convenienc^ 



B. Interference metrics and scheduling assumptions 

The interference between IS and AH networks is modeled 
by the average temporal overlap between both systems. Based 
on the sensing result at the beginning of the frame and knowl- 
edge of the CTMC parameters of the AH links, transmission 
time and power are allocated. 

The allocation of transmission time consists of specifying 
duration and placement of the transmission within the current 
frame. We first show that it is optimal to transmit at the 
beginning (the end) of the frame if the sensing outcome is 
idle (busy). 

Lemma 1: Assume that a p fraction of transmission time 
needs to be allocated to a sub-channel, on which the AH 
user's ON/OFF behavior is modeled by the CTMC ©■ Based 
on a sensing outcome at the beginning of the frame, the 
minimum expected overlap with the ON period of the AH 
user is achieved by 

• transmitting at the beginning of the frame {i.e., during 
[0, pT]) if the sensing outcome was idle and 

• transmitting at the end of the frame {i.e., during [(1 — 
p)T, T] if the sensing outcome was busy. 

Proof: see appendix. 
Based on Lemma 1, we derive the expected time overlap 
between IS and AH transmissions, conditioned on the sensing 
result y e {0, 1} at the beginning of the frame. Consider a sub- 
channel 71, which overlaps with AH band i = g{n). Then, the 
activity of AH user i is given by the CTMC > 0} 

with parameters A; and pi. Transmitting for a p fraction of 
the frame, leads to the expected time overlap 



0n,o(p) 



|E|y'l{x^(^)=ijdeX(0) = 

pT 

= ^J Pr(X(C)==l|X(0) = 0)de (3) 



if the sensing result was idle. By substituting ([T]i it is then 
easy to show that 

0«.o(p) = - 1)) . (4) 

In the case of a busy sensing result we obtain 



(P) 



(5) 



This derivation made use of the fact that the sensing results 
of two IS sub-channels are either perfectly correlated (if they 
overlap with the same AH band) or statistically independent 
(if they overlap with different bands). Therefore, the prediction 

^The above formulation encompasses a channel capacity formulation (for 
K = 1) as well as the case of variable-rate M-QAM in which case k = 
1.5/(— In BER) is chosen such that a target BER is met 1211 . 
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performance of a specific sub-channel cannot be improved by 
using sensing results from other sub-channels. 

Lemma 2: The functions (j>„.o{p) and (j)n,iip) are strictly 
convex and increasing in p. 

Proof: Both (f>n.o{p) and 0„,i(/o) are nonnegative linear 
combinations of a convex and a strictly convex function. 
One is linear, the other an exponential function with nonzero 
exponent. The monotonicity can easily be verified by differ- 
entiation. 

III. Optimal Frame-Level Allocation 

Consider a single IS client, which minimizes the time over- 
lap between IS and AH transmissions subject to maintaining 
a rate constraint across the IS channel to the IS base station. 
Mathematically, this leads to problem PI 



min y^(/>n,y„(Pn) 
p p '-^ 



S.t. ^p„log(l + ^) >i? 



Pn 



Pn<P 



(6) 



(7) 



(8) 



p„ > 0, 1 < n < A/- (9) 
< /0,i < 1, l<n<N, (10) 

with rate constraint (|7]i and power constraint ([8]l. It is straight- 
forward to show that PI is a convex optimization problem 
since the objective function is convex (by Lemma 2), the rate 
constraint (once rewritten in standard form) is convex by the 
perspective property [22], and all other constraints are linear. 

A solution to PI can be found by general solution tech- 
niques in polynomial time 1221 . For this specific problem, 
however, it is possible to show a special structure that enables 
us to gain further insight into the problem. 

A. Optimality conditions and solution structure 

The solution structure is obtained by introducing Lagrange 
multipliers 7 and e for the rate and power constraint, respec- 
tively. This leads to the Lagrangian 



n 

\ Pn J 



+ e 



P 



(11) 



The Karush-Kuhn-Tucker (KKT) optimality conditions are 
then given by the constraints (fTli-dToli of PI, nonnegativity 
constraints for the Lagrange multipliers, 7 > 0, e > 0, the 
slackness conditions 



7 



i?-^p:iog(i 



= 



the condition 

9i(p,/3;7,e) 



d-pn 



= 0, K>o 
>o, K = o 



(12) 
(13) 

(14) 



i9i(p,p;7,e) 



dp-n 



= 

€(0,1) 
= 1 



(15) 



and 

> 0, pI 
= 0, pi 
< 0, p* 

Expressions fT4i and ( fTSl l can be understood on an intuitive 
level by noting that for p* and p* to be minimizers of 
L{p, p; 7, e), their partial derivative must equal zero unless 
they lie on the boundary of the feasible set. 

Solution structure for p*^: By substituting (fTTT l into (fT4l i and 
solving for p* we arrive at 



(16) 



where ly := 7/e has been introduced to simplify notation in 
what follows. For any fixed value of pn, (O represents a 
water filling solution 1231 . 

Solution structure for p* .• The optimal transmission time 
allocation is obtained by substituting ( fTTT l and ( fT6l ) into ( fTSl ). 
For an idle sensing result, ?;„ = we obtain. 



Pn = 



(Ai+Mi)T 
1, 



log- 



(17) 

where Coa = Xi/{Xi + pi){l-exp{-{Xi + pi)T) and i = g{n) 
denotes the AH sub-band that overlaps with sub-channel n. In 
the above equation we have defined 

{^Pn " 1) + 



hn{v) [log(i//3„)]" 



(18) 



1 + {ly/Sn - 1)+ 
to simplify notation. Similarly, we can obtain the solution 
structure for the case of a busy sensing result, i/„ = 1, 



0, 

u 
1, 



log 



Cl,t < lhn{v) < 1 

jhn{iy) > 1, 



(19) 

where Ci.i = Xi/{Xt+ Pi)(l + Pt/\iexp{-{X^ + pi)T). Note 
that these closed-form expressions depend on the Lagrange 
multipliers only through the term ^hn{v) which does not 
depend on the AH activity parameters A; and pi. Further, 
the transmission time allocations (fT7])-(fT9]) are monotonic with 
respect to this term. 

B. Iterative solution algorithm for 7 and v 

To find the optimal power and transmission time allocation 
based on the above closed-form expressions, we present an 
algorithm for finding the pair [7*,;/*], which corresponds to 
the optimal solution of PI. 

For any pair [7, v\ the power allocation p(7, v) and trans- 
mission time allocation p(7, v) define the optimal solution to 
PI with modified rate constraint 



n 

and modified power constraint 

-P(7, V) := ^ p„(7, v)[v - ^ 



(20) 



(21) 
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where p„(7, j/) is given by ( fTTT i or ( fT9] l (depending on the 
sensing resuh). The fact that this solution is optimal for 
rate constraint i?(7, ly) and power constraint P{j, v) follows 
directly from the KKT optimality conditions, which are nec- 
essary and sufficient for convex optimization problems 1221 . 

Based on the above, finding the pair [7*, v*] corresponding 
to the given rate constraint R and power constraint P could 
theoretically be performed by searching all pairs [7, i^]. In this 
section we show, however, that i?(7, v) and P(7, v) exhibit 
some monotonicity which enables us to use the bisection 
method for finding [7*,;/*] with guaranteed convergence and 
low complexity. We first study the case of keeping 7 fixed and 
adjusting v such that the rate constraint is met with equality. 
Then, we show that the allocated sum power decreases with 
7- 

Adjusting v to meet the rate constraint: We first consider 
the case of adjusting v such that -R(7, v) = R while keeping 7 
fixed. It is easy to verify that for any n, hn{v) is nondecreasing 
in V. Therefore, /9„(7,i/) increases with 7 as well as can be 
seen from ( [TtI i and ( fT9] l. Further, since /0„(7, v) increases with 
V for fixed 7, so does i?(7, v). We can exploit this property 
to find the v for which i?(7, v) = Rhy the bisection method. 
First, we can find upper and lower bounds, v.^ and vi, for this 
value. These bounds are guaranteed to exist since i?(7, v) 
00 for ly 00 and R{j, i^) for ly 0. Once these bounds 
have been obtained the bisection method iteratively finds i/* 
with guaranteed convergence. 

Adjusting 7 to meet the power constraint: Having obtained 
an algorithm for finding v for arbitrary 7 such that the rate 
constraint is satisfied, we study the behavior of the power 
constraint as 7 is adjusted. As 7 is varied, we continue 
to adjust V such that the rate constraint is satisfied at all 
times. The pair of Lagrange multipliers is therefore given by 
[7,'^*(7)]- 

The slackness conditions imply that at the optimal solution 
both rate and power constraints are met with equality. From 
(|20] | we observe that decreasing 7 requires increasing v in or- 
der to continue meeting the rate constraint. Further, decreasing 
7 reduces the objective function because \og{vj3n) > log(j^/3ri) 
for V >v, enables us to reduce Pn{l, for at least some 

n. 

Since decreasing 7 requires increasing v, the allocated sum 
power increases as 7 decreases. Intuitively, a constant rate 
constraint ( |2Q] | requires ( 1211 1 to increase because the term 
[zy*(7) - increases faster than [log(z^*(7)/3„)]+. The 

following lemma shows this rigorously. 

Lemma 3: The sum power J^nPn associated with alloca- 
tion [7, J'* (7)] is a decreasing function in 7. 
Proof: see appendix. 

Lemma [3] enables us to find 7* again by the bisection 
method. Assuming that PI is feasible which we will as- 
sume hereafter, there exist bounds ju and 7; such that 
P{lu,v*hu)) < P < P{-iuv*{ii)). Therefore, by starting 
the bisection method from these points we can find the pair 
[7*,z/*] with guaranteed convergence. The solution algorithm 
is shown in detail in Fig. |2] The inner loop (lines 4-14) 
correspond to finding v*{'^), whereas the outer loop finds 7*. 



Algorithm 1: Solution Algorithm 



1 Initialization. Obtain bounds vi, v^, 7;, 7„; 

2 repeat 

3 7 ^ (7« - 7;)/2; 

4 repeat 

5 i> ^ - f/)/2; 

6 Find time allocation p„(z>) using (flTl l; 

7 Find power allocation Pn{j, v) using ( fTSI l; 
Compute achievable rate r(p„,p„) using (|2]i; 
if r{pn,Pn) > R then 

I i> 
else 



9 

10 
11 
12 

13 
14 
15 
16 
17 
18 
19 
20 
21 



I n ^ V 
end 

until < i? - r{pn, p„) < ; 
Find time allocation p„(7, i>) using (flTl l: 
Find power allocation Pnil, i^) using ( fT6b : 
if J2nPn then 
I 7i ^ 7 
else 

I 7« ^ 7 
end 



22 until < P - Y^nPn < ep 



Fig. 2: Algorithm for finding the optimal Lagrange multipliers 
7 and ly for problem (l6]l-(fT0li. The inner loop (lines 4-14) find 
J/* (7) which satisfies the rate constraint (|7]i. The outer loop 
determines 7*, which satisfies the power constraint dHJ. 



C. Properties of optimal allocations 

Beyond simplifying solution algorithms, the structured solu- 
tions also enable us to make some qualitative statements about 
the optimal resource allocation. 

Ordering in terms of IS channel quality: We first investigate 
the ordering with respect to the IS channel coefficients /3„. 
This scenario is shown in Fig. [3a] for iV = 4 sub-channels 
with coefficients (3 — [.9, 1.1, .5, 1.5]-^. We observe that, for 
any rate constraint, more transmission time is allocated to the 
sub-channel with higher channel coefficient, i.e., /3j > jSj => 
Pi > Pj provided all other parameters and the sensing results 
are identical. On an intuitive level, this result captures the fact 
that in channels with high Pi we can achieve the same rate in a 
shorter transmission duration using the same amount of power 
Mathematically, the result follows from the monotonicity of 
the optimal solution and hn{v). 

Ordering in terms of sensing results: Similar to the IS 
channel, the optimal transmission time allocation can be 
ordered with respect to the sensing results. If channels have 
the same IS channel coefficient, (3i = f3j but different sensing 
results then it is preferable to allocate more transmission time 
to the idle channel, i.e., pi > pj. This is illustrated in Fig. [3b] 
for AT = 4, /3 = [.9, .9, 1.1, 1.1]^, and y = [0,1,0,1]^. It is 
also interesting to note that some transmission time is allocated 
to frames with busy sensing results even when idle frames are 
not yet used to the maximum extent. 
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rate constraint R [bps] rate constraint R [bps] 

(a) Ordering by IS channel quality (b) Ordering by sensing outcome 

Fig. 3: Structure of the optimal transmission time allocation for varying rate constraint. The solution can be ordered with 
respect to IS channel coefficients or sensing outcomes. 



IV. Optimal Average Resource Allocation 

Problem formulation PI required that rate and power con- 
straint are met in every frame, even if sensing outcome or 
IS channel quality are disadvantageous. In practical systems, 
satisfying rate constraints at the frame-level is usually unnec- 
essary; it suffices to maintain average rate constraints across 
time. This less stringent requirement can be used to further 
reduce interference by allocating less transmission time during 
frames with adversarial channel/interference conditions, while 
compensating for the rate decrease during frames with better 
conditions. Ultimately, this leads to an improved resource 
allocation across both frequency (the sub-channels of the IS 
system) and time (consecutive frames of the IS system). 

This section introduces such an average rate formulation by 
averaging across the temporal activity of the AH network and 
random IS channel coefficients. Further, this section introduces 
two reference schemes that help to put the performance of the 
optimal resource allocation in perspective. 



yG¥ n 

Pn,y > 0, Vy e Y, 1 < 77, < TV (25) 
< pn.y < 1, Vy e Y, 1 < 7i < iV, (26) 

where -q^^o = fJ.i/{K + fJ-i), = + f^t), and due 

to the independence of the AH sub-bands, rjy = Ilf^i Vi,yi- 
Note that this optimization problem has 2*^ as many decision 
variables because power and transmission time allocation may 
be different for every possible sensing outcome. The fact that 
the decision variables grow exponentially with M is not of 
major concern, because M (the number of parallel AH bands) 
is typically quite small (in the order of one to five). 

Problem P2 can be solved similar to problem PI. In 
particular, by forming the Lagrangian, introducing Lagrange 
multipliers 7 and e, and taking the derivative with respect to 
the decision variables, we obtain a similar solution structure 
as in the frame-level problem. 



A. Formulation and solution structure 

The average rate formulation requires associating probabili- 
ties with all possible sensing outcomes. While there are a total 
of N sub-channels available, the sensing outcomes for sub- 
channels that overlap with the same AH band will be identical. 
Therefore, for M sub-bands, there are a total of 2*^ possible 
sensing outcomes. Let the set of all possible sensing outcomes 
be represented by Y = {0, 1}^^ where y [yi, . . . , t/a/]^ G Y 
denotes the sensing outcome per sub-band. 

Problem P2 of optimally allocating power and transmission 
time then becomes 

min 5I'^yH'^">!^s(")(^">y) (22) 

Pn'.y yGY n 

s.t. Y.^yll P-.y log f 1 + ^ ^ (23) 



B. Reference schemes 

Having incorporated random sensing outcomes into our 
formulation, we introduce two reference schemes in order 
to put the performance of the optimal resource allocation in 
perspective. 

No sensing case: As a first benchmark, consider an ap- 
proach that allocates power but does not perform any transmis- 
sion time optimization. This case corresponds to conventional 
resource management in IS systems, which simply allocates 
power to make the best use of the uplink channel. Mathe- 
matically, this is formulated as minimizing pn subject to 
the constraints (l7]i-(fT0]i. We assume that for any sub-channel 
with Pn > 0, the sub-channel is used for the entire frame 
duration. Unused carriers for which p„ — are not allocated 
any transmission time. 
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Idle-frame allocation: Another possible reference scheme 
performs spectrum sensing but allocates resources in a sub- 
optimal way. Specifically, consider allocating the entire frame 
by setting /9„ = 1 for all idle sub-channels while completely 
avoiding busy sub-channels by setting p„ = for all n with 
Un = 1. In the average rate formulation, this method can be 
formulated mathematically as minimizing J2yeY^y^nPn-y 
subject to (l23Tl-(l26]l and the additional constraint that sub- 
channels with busy sensing outcome are never allocated. Note 
that the above optimization problem may be infeasible even 
when P2 is feasible because we are imposing the additional 
restriction of never transmitting during busy frames. To ensure 
that the reference scheme is always feasible when P2 is, we 
force allocation to busy channels if the resulting optimization 
would otherwise be infeasible. 



C. Allocation for random IS channels 

The previous section considered average rate constraints 
with respect to temporal AH activity but fixed IS channel 
/3. This section further extends the analysis to the case of 
random IS channel coefficients (3. The optimization problem 
P3 is given by 

f% / E E (/^"(y' f^))dF{f3) (27) 

s-t- / E'?yE/^«(y'/3)iog(i + 2^^||^)dF(/3)>i? 

(28) 

/ E^yE^^»(y'^)^^(/3) (29) 



p yGY n 

Pn(y,/3)>0, Wy eY,l<n<N (30) 
0<p„(y,;3)<l, VyeY,l<n<iV,, (31) 



where the decision variables p(y, f3) and p(y, /3) correspond 
to the power and transmission time allocation that is used for 
sensing outcome y and IS channel condition f3 and F{/3) is 
the cumulative distribution function of /3. By again forming 
the Lagrangian and computing the derivative with respect to 
the decision variables, it is easy to show that the structured 
solutions ( fT4b and (fT7]i-(fT9]l again hold. Therefore, it is again 
possible to express the allocation as a function of the Lagrange 
multipliers [7, ly]. The rate constraint can then be evaluated by 

/ E'?yE/'"(y"3)[i°g('^/3")]+rf^(/3) (32) 

and the allocated sum power is given by 

/E'?yE'^"(y'/3)(^-i)^'^^(^)- (33) 

p ySY n 

While the above integrals can only be evaluated numerically, it 
is possible to again find the optimal solution via the bisection 
method. 

D. Numerical results 

This section presents numerical performance results for the 
optimal average rate resource allocation and compares them to 
the reference schemes introduced in this section. The results 
were obtained for iV = 5 sub-channels and a single AH sub- 
band M ~ 1. The prediction parameters were A = /i = ls^^ 
and the IS channel coefficients were flat Rayleigh fading and 
statistically independent. We further assume a block fading 
scenario in which the IS channel varies slowly compared to 
the frame duration. 

The performance for fixed IS channel and random sensing 
results is shown in Fig. |4a| for T = 1 s and in Fig. |4b] for 
T = .Is. The plot shows the average transmission overlap 
between IS and AH network versus the achieved IS rate 
(note that the achieved IS rate and not the rate constraint 



SUBMITTED TO IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS 



8 



is plotted). The performance results are averaged over 100 
realizations of the IS channel. Since the IS channel is not 
modeled statistically, it is inevitable that for some realizations 
of f3 problem P2 is infeasible. The outage probability, which 
is identical for all three schemes, is therefore shown in Fig.l4cl 
to put the results in perspective. Typical outage probabilities 
of approximately 10% correspond to an IS rate of about 
0.7 bps. At this rate, plots (a) and (b) show that a significant 
performance gain is achieved by performing sensing-based 
transmission time allocation. 

The performance ordering reflects our expectations. The 
idle-frame allocation scheme outperforms the no-sensing case 
but shows a quite significant performance gap with respect to 
the optimal allocation, especially for low IS rates. Further, all 
curves show increasing interference as the IS rate increases. 
This is expected, since high IS rates prevent the IS system from 
being able to accommodate the AH links. The plots also show 
that idle-frame allocation and no-sensing scheme converge for 
high IS rates, because allocating only idle frames is almost 
always infeasible (and therefore busy frames typically need to 
be used as well). 

By comparing Fig. ga] (T = 1) and Fig. gb] (T .1), we 
observe that while the performance of the optimal scheme does 
not change significantly, the idle-frame reference performs 
much better This is intuitive, because by reducing the frame 
length, it is easier to "fill up" the idle periods of the AH 
network. The performance of the no-sensing scheme remains 
unaltered and is the same in both figures. 

The performance for average IS channel coefficients and 
random AH behavior is shown in Fig. |5] which compares 
the solution of P3 with the same reference schemes. We 
can observe that by exploiting the channel variability and 
allocating across frequency and time, we can further reduce 
interference. Otherwise, the performance trends are similar 
to those of Fig. |4] Note that idle-frame allocation does not 
achieve the same channel capacity as the optimal scheme 
because it only transmits in frames with an idle sensing result. 

V. Allocation for Multiple IS Users 

The previous section derived the optimal power and trans- 
mission time allocation assuming that an orthogonal set of 
sub-channels had already been assigned to each IS user. This 
enabled us to consider each of the terminals individually and 
perform resource allocation based on local sensing results. 

In practice, the IS base station needs to assign sub-channels 
to each of the IS users without knowing what the sensing 
outcomes will be. We therefore consider the problem of 
optimal sub-channel allocation based on minimizing average 
interference metrics. This leads to a similar formulation as 
compared to the average interference case in Sec. |IV] Once 
a sub-channel allocation has been computed and fed back 
to the IS users, they can use the locally available sensing 
results to optimize their medium access. The optimal sub- 
channel allocation is a combinatorial problem, which is com- 
putationally more challenging than problems PI through P3, 
which could be analyzed based on convex optimization. While 
a general analysis of this problem goes beyond the scope 
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Fig. 5; Performance result for random IS channel coefficients. 



of this paper we find the optimal solution for fairly small 
problem instances by exhaustive search. A comparison with 
heuristic allocations suggests that efficient greedy sub-channel 
allocation algorithms developed for related problem setups can 
be adopted to this problem and yield a performance close to 
optimal. 

A. Optimal sub-channel allocation 

The problem of optimal sub-channel allocation involves 
assigning orthogonal sets of sub-channels to each terminal, 
such that the overall interference is minimized; see Fig. |6] 
A mathematical formulation can be based on problem P2. 
Specifically, define /(A) as the optimal value of P2 where 
the summations over sub-channels are restricted to n G A,;.e., 
the sub-channels on which a specific IS user operates. Define 
/(A) := oo if P2 is infeasible. 

Assume that the IS base station is serving a total of U 
users. The problem of optimal sub-channel allocation is then 
formulated as 



min V/(A„) 
{A„} ^ 



u=l 



s.t. AiU---UA[/ = {!,..., TV} 



n. 



(34) 

(35) 
(36) 



where A^ denotes the set of sub-channels assigned to termi- 
nal i. Due to i35[ and ( |36] | the sub-channel allocations are 
mutually exclusive and collectively exhaustive. Note that the 
above problem does not require knowledge of the sensing 
outcome at the individual terminals. The base station only 
requires knowledge of the CTMC parameters A and /i, as well 
as knowledge of the IS channel coefficients (3. 

The above problem is difficult to solve due to its combina- 
torial nature and conventional sub-channel allocation methods 
are not easily extended to incorporate the additional dimension 
of allocating transmission time. For small problem instances. 
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however, the optimal allocation can be found by exhaustive 
search. 



B. Suboptimal algorithm 

The problem of optimal sub-channel allocation in multicar- 
rier systems has been well-studied in the absence of trans- 
mission time allocation. Standard methods typically minimize 
the total transmit power subject to rate constraints. In our 
setup, this can be formulated mathematically as minimizing 
Sm=i Sn6A subjcct to rate and power constraints for the 
individual terminals. While the resulting optimization problem 
is still combinatorial, efficient approximation techniques have 
been developed with close-to-optimal performance. 

Allocating sub-channels in this manner, can be used as an 
effective heuristic. Since good channel quality results in lower 
average transmission time, we conjecture that conventional 
sub-channel allocation may be a good approximation to the 
optimal interference-aware sub-channel allocation. Numerical 
results show that this is indeed the case in the scenarios we 
have examined. 

C. Numerical results 

Numerical results for the multi-terminal case are shown in 
Fig. |2] The total average transmission time overlap (summed 
over all IS users) is plotted with respect to the rate constraint 
for each individual IS users (constraints are assumed to be 
identical). The performance trends are the same as in the 
case of a single IS user For low rate requirements we 
can effectively mitigate interference by assigning resources 
judiciously. On the other hand, as rate requirements become 
more stringent, there is less flexibility in accommodating the 
AH links. The scenario plotted in Fig. [7] corresponds to U ~ 3 
terminals, N ~ 5 sub-channels, and flat Rayleigh fading IS 
channel coefficients. 

The performance of the optimal and suboptimal sub-channel 
allocation schemes is very similar regardless of the rate 
constraint. This suggests that minimizing the total transmission 
power is a reasonable approximation to the optimal sub- 
channel allocation. In future work, we plan to corroborate this 
conjecture in more detail. 

VI. Conclusion 

In conclusion, this paper has introduced a novel cognitive 
coexistence framework, which enables infrastructure systems 
to coexist with local, ad-hoc or peer-to-peer communication 
links. Based on sensing and predicting the interference patterns 
of these ad-hoc links, the infrastructure system allocates power 
and transmission time judiciously such that overlaps with the 
ad-hoc links are minimized. 

We analyzed the problem within the framework of convex 
optimization and derived closed-form solutions at the frame- 
level. These results were extended to the average rate case to 
reduce interference further by allocation across both frequency 
and time. Finally, we addressed the case of multiple IS users 
and provided more insight on how sub-channel allocation can 
be performed. 
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cluster 1 cluster 2 cluster U 

Fig. 6: System setup for the multi-terminal case. 
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Appendix A 
Proof of Lemma I 

We denote the IS transmissions within the current frame by 
a finite set of closed and disjoint time intervals = [a^j^A;] 
where each Ife C [0, T\ corresponds to a contiguous trans- 
mission of the IS user. Clearly, this formulation incorporates 
possible pauses between IS transmissions. We also require 
^Yuh bk ~ dk = pT, because a total of pT transmission time 
needs to be allocated. 

First, consider the case of an idle sensing result at the 
beginning of the frame, say at time < = 0. Then according 
to ^ and (HJ, the expected time overlap is given by 



Pr(X(0 = 1|X(0) = 0)d^ 



1 

T 



E 



d^. (37) 
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Since the integrand is strictly increasing in p, the above 
expression is minimized by transmitting contiguously during 
the time interval [0,pr]. 

In the case of a busy sensing result, an equivalent approach 
leads to a strictly decreasing integrand and therefore it is 
optimal to transmit during the time interval [{\ — p)T,T] in 
that case. 

Appendix B 
Proof of Lemma 3 

In Sec. lIII-Bj we have defined the sum power corresponding 
to the pair of Lagrange multipliers [7, v] as i-'(7, v). Further, 
we showed that by keeping 7 fixed and varying v it is possible 
to find a V* (7) for which the rate constraint is satisfied with 
equality. To simplify notation let us now define P(7) as the 
sum power associated with [7, v*{'-^)]. 

The proof that P(7) decreases with 7 proceeds by contra- 
diction. First, we note that 7^0 implies P(7) 00 due 
to the structure of the optimal solutions (fT7li-(fT9]). Assume 
now that P(7) is not monotonically decreasing. Then, because 
P(7) is continuous, there exist two different values of 7, say 
71 and 72, such that P(7i) = ^(72)- 

Based on the KKT conditions stated in Sec. IIII-AI it is 
easy to verify that both 71 and 72 correspond to optimal 
solutions of Problem PI with rate constraint i?(7i, i^*(7i)) and 
power constraint P(7i, i/*(7i)). Further, from the structure 
of the optimal solutions it is clear that the transmission time 
allocations associated with 71 and 72 must be different, that 
is, p(7i) ^ p{l2)- This is a contradiction, however, because 
Problem PI has a strictly convex objective function and 
therefore at most one optimal solution. 
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